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ABSTRACT A bioinformatic and phylogenetic study has been performed on a family of penicillin-binding proteins 
including D-aminopeptidases, D-amino acid amidases, DD-carboxypeptidases, and p-lactamases. Significant homol- 
ogy between D-aminopeptidase from Ochrobactrum anthropi and other members of the family has been shown and 
a number of conserved residues identified as S62, K65, Y153, N155, H287, and G289. Three of those (Ser62, Lys65, and 
Tyrl53) form a catalytic triangle - the proton relay system that activates the generalized nucleophile in the course of 
catalysis. Molecular modeling has indicated the conserved residue Lys65 to have an unusually low pKa value, which 
has been confirmed experimentally by a study of the pH-profile of D-aminopeptidase catalytic activity. The resulting 
data have been used to elucidate the role of Lys65 in the catalytic mechanism of D-aminopeptidase as a general base 
for proton transfer from catalytic Ser62 to Tyrl53, and vice versa, during the formation and hydrolysis of the acyl- 
enzyme intermediate. 
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INTRODUCTION 

D-aminopeptidase from Ochrobactrum anthropi possesses 
a unique structural organization, high stereospecificity, and 
shows catalytic activity towards a wide range of D-alanine de- 
rivatives [1]. It is a member of the serine hydrolases superfami- 
ly and acts as a homodimer, with each subunit consisting of 
three structural domains. One of those - a catalytic domain - 
has the so called |3-lactamase fold [2]. Sequence comparison has 
revealed a strong evolutionary relationship of D-aminopepti- 
dase with DD-carboxypeptidase and |3-lactamase [3]. 

The catalytic mechanism of D-aminopeptidase has not 
been discussed in the literature; however, suggestions con- 
cerning other enzymes of the family - D-amino acid amidase 
from Ochrobactrum anthropi, as well as -lactamases and pen- 
icillin-binding proteins, can be taken into account. Different 
views on the catalytic mechanism of penicillin-binding pro- 
teins are presented, and several residues in close proximity 
to the catalytic serine are considered as potential candidates 
for the role of general base in the course of the enzymatic 
reaction [4- 6]. While choosing between them, it is important 
to look at the pH-profile of enzyme activity and possible pK 
shifts of the residue due to its environment, since the general 
base is bound to act in deprotonated form. 

It is important that molecular modeling of the deacylation 
step has been performed for reactions catalyzed by DD-pep- 
tidase R61 from Streptomyces sp. and class C |3-lactamase P99 



from Enterobacter cloacae by means of QM/MM methods [7]. 
It revealed the leading role of the active site tyrosine residue 
in the deacylation of the catalytic serine of those enzymes. 

Yet, an analysis of the literature shows that a comprehen- 
sive view on the catalytic mechanism of D-aminopeptidase is 
still lacking. In this work, we aim to use bioinformatics and 
molecular modeling to elucidate the role of the Lys65 residue 
in the catalytic triad of D-aminopeptidase from Ochrobac- 
trum anthropi. 

MATERIALS AND METHODS 

Experimental study 

D-aminopeptidase was purified according to a procedure 
described earlier [1]. The colorimetric substrate - D-alanine 
p-nitroanilide (D-Ala-pNA) - was produced by Bachem. 
Tris(hydroxymethyl)methylamine (Tris) of research-grade 
purity produced by SERVA Electrophoresis was used to pre- 
pare buffer solutions. 

Kinetic assays were performed using progress curve anal- 
ysis. The studies were carried out in a spectrophotometric 
cuvette of 500 ul and optical path length of 1 cm thermostat- 
ed at 25°C. Enzymatic hydrolysis was initiated by adding a 
small amount of the enzyme solution to the reaction mixture 
containing the substrate in a concentration approximately 4 
times higher than its K M . Changes in absorption were reg- 
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Fig. 1 . Phyloge- 
netic tree based on 
structural alignment 
of penicillin-binding 
proteins and 
primary structure 
motives containing 
important catalytic 
residues: Ser62 
and Lys65 (SXXK), 
Tyr153([YS]XN) 
and His287 ([KH] 
XG). Leaves are 
named according 
to PDB databank 
accession numbers. 
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istered by a Shimadzu UV-1601 spectrophotometer in the 
Kinetics mode at 450 nm. The final levels of absorption were 
not higher than 2 in all of the assays. 

To maintain the set pH value in the reaction mixture, the 
enzymatic reactions were carried out in a 0.1 M Tris-HCl 
buffer. A Hamilton Slimtrode pH-sensitive electrode was 
used in the preparation of the buffer solutions. 

The progress curves obtained were processed by data lin- 
earization in t/lnfp^/fp^-p)) - p/lnfp^/fp^-p)) coordinates, 
where t is time, p - current concentration of the product, and 
p^ - the final concentration of the product. This anamorphosis 
allows to determine the value of the K../V and 1/V ratio 

M' max ' max 

as the line's slope and y-intercept, respectively. The nonlinear 
regression of the V max /K M dependence on pH and the com- 
putation of experimental pKa values were performed using 
SciDAVis software [8]. 

Bioinformatics 

Homology search. In all homology search procedures, a se- 
quence or structure of D-aminopeptidase from Ochrobactrum 
anthropi (PDB entry 1EI5) was used as a query. 

The sequence-based homology search was carried out 
using the PSI-BLAST [9] algorithm v. 2.2.18 to scan a "non- 
redundant" protein sequence database. The resulting sample 
was filtered with a 95% pairwise identity threshold to elimi- 
nate redundancy and then aligned using t_coffee [10], mafft 
[11], and probcons [12] joined together by the consistency- 
based statistics implemented in t_coff ee. 

The structure-based homology search was carried out by 
scanning the PDB protein structure databank using the SSM 
[13] procedure. Hits were discriminated in case of high sec- 
ondary structure elements mismatch with the 1EI5 structure. 
The resulting sample of three-dimensional structures was 
aligned using the MUSTANG [14] software. 



Bioinformatic analysis. The phylogenetic analysis of se- 
quence and structural alignment was performed using the 
phylip package [15]. The phylograms were constructed using 
distance-based methods with the neighbor-joining algorithm. 
The bioinformatic analysis was carried out using the original 
ZEBRA v. 3.2 software with a statistical threshold level of 
2.2 xlO" 43 . 

Visualisation. The Jalview [16] program was used to look 
through multiple sequence alignments. Visualization of three- 
dimensional structures and a structure-based multiple align- 
ment was done using PyMol [17]. Generation of phylogenetic 
trees was done using phylip [15]. Generation of sequence pat- 
terns logotypes was done with the WebLogo [18] Internet 
service. 

RESULTS AND DISCUSSIONS 

Bioinformatic analysis of enzymes homologous 
to D-aminopeptidase 

Data related to the penicillin-binding protein family includ- 
ing D-aminopeptidases, D-amino acid amidases, and alkaline 
D-peptidases were collected and analyzed. The UniProt pro- 
tein sequence database and PDB protein structures databank 
were screened a priori to identify all significant sequences and 
structure-based homologs of D-aminopeptidase. The result- 
ing set of 734 sequence homologs and 24 structural homologs 
was sampled and filtered to acquire the most informative set. 
As a result of structural alignment, a significant similarity of 
the active site regions of D-aminopeptidase, alkaline D-pepti- 
dases, D-amino acid amidases, and |3-lactamases was shown 
for both the sequence and structure levels (fig. 1). 

The following residues were identified as conserved in the 
active site of penicillin-binding proteins using the original 
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Fig. 3. Stereo view of the D-aminopeptidase active site and 
spatial organization of catalytic triad. The S62, K65, Y153, N155, 
H287, and D225 residues are highlighted. 



Fig. 2. An active site fragment of the structural alignment of 
penicillin-binding proteins - structural homologs of D-ami- 
nopeptidase from Ochrobactrum anthropi. Location of four con- 
served active site residues — Ser62, Lys65, Tyr153, His287 — 
conserved in D-aminopeptidase, as well as in other members of 
the family, is shown. 



ZEBRA software developed in our laboratory: 287H, 153Y, 
155N, 289G, 273G, 293G, 270Y, 65K, 224G, 62S, 68T, 294W, 
64S, 151Y, 2281, 601, and 288G (numbered according to the 
1EI5 structure and sorted in decreasing significance). Con- 
sidering that the conservation of a residue in a protein struc- 
ture indicates an evolutionary pressure on that position and 
thus underlines its functional or structural importance [19], 
we suggest that such residues are important to the D-ami- 
nopeptidase catalytic mechanism. A common alpha-beta do- 
main identified as a "three-layer sandwich" by CATH [20] 
structure classification was shown to contain active site resi- 
dues in all studied penicillin-binding proteins (fig. 2). 
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D-aminopeptidase structure analysis 

The D-aminopeptidase structure (PDB entry 1EI5) analysis 
reveals a pair of amino acid residues, Tyrl53 and Lys65, lo- 
cated in near proximity and in approximately equal distance 
to the catalytic Ser62's Oy atom. All three residues form a 
nearly equilateral triangle in the enzyme active site (fig. 3). 
Such a location of residues implements a special organiza- 
tion of the proton relay system when the hydrogen atom of 
the serine's hydroxyl-group is directed toward the center 
of the triangle and shared among all of the residues of the 
catalytic triad. The specific organization of this catalytic tri- 
angle is based on the irregular properties of the Lys65 residue 
capable of accepting a proton in neutral and slightly alkaline 
media, which lends high reactivity to the Ser62 residue at the 
formation of the acyl-enzyme intermediate. 

A considerably lower pKa value of the Lys65 residue equal 
to 7.8, compared to the ionization of regular lysine residues 
in proteins with pKa 10-11, was observed upon calculation 



Fig. 4. pH-profile of D-aminopeptidase catalytic activity. The 
theoretical curve was calculated according to the equation v = 
v 0 /(1 +[H + ]/K 1 +K 2 /[H + ]), where v Qf pK 1 and pK 2 values are 
equal to 0.53, 7.4 and 8.8, respectively. 



of the ionization properties of D-aminopeptidase active site 
residues by the PROPKA QSAR method [21]. A high-eval- 
uated pKa value of the Tyrl53 residue equal to 11.85 should 
be noted as well. 

Experimental pH-profile of D-aminopeptidase catalytic 
activity 

Experimental data are in good agreement with the results 
obtained by molecular modeling and shed light on the role 
of the Lys65 residue in the functioning of the catalytic triad. 
The pH-prof ile of the D-aminopeptidase catalytic activity has 
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Fig. 5. Schematic presentation of acylation (a) and deacylation (b) of the Ser62 residue in the catalytic mechanism of D-aminopepti- 
dase: (a) organization of catalytic triad and its role in the formation of the first tetrahedral intermediate followed by the formation of 
acyl-enzyme are shown; (b) nucleophile (water molecule) binding, activation, nucleophilic attack followed by formation of the second 
tetrahedral intermediate and regeneration of a free enzyme are shown. 



a bell-shaped form with a pKa value of 7.4 and a pKb value 
of 8.8 (fig. 4). On the basis of the molecular modeling and the 
experimental study, the amino acid residue with a pKa of 7.4 
was referred to Lys65. 

Proposed catalytic mechanism of D-aminopeptidase 

Just as in the case of other serine hydrolases, the reaction cata- 
lyzed by D-aminopeptidase follows a 3-step kinetic scheme 
with the formation of a covalent acyl-enzyme intermediate and 
its subsequent hydrolysis or transfer of the acyl group to an 
external nucleophile. Because of the extremely low pKa value 
of its terminal amino group, the Lys65 residue plays a specific 
role in the D-aminopeptidase catalytic mechanism (fig. 5): 

Being uncharged at the pH-optimum of the enzymatic 
reaction, the Lys65 residue acts as a general base, while the 
Oy atom of Ser62 attacks the carbonyl group of a substrate: 
Lys65 assumes a proton from the attacking OH-group when 
the first tetrahedral intermediate is formed at the acylation 
step. 



At the decomposition of the first tetrahedral intermediate 
followed by the formation of the acyl-enzyme and release of 
the first reaction product, its leaving group gathers a pro- 
ton donated by the OH-group of Tyrl53, whose acidity is in- 
creased due to the proximity of a positively charged Lys65 
residue, and at the same time the formed oxyanion of Tyrl53, 
being stronger base, captures a proton from the Lys65 resi- 
due. 

At the deacylation step, the water molecule (or molecule 
of another nucleophile) is activated through the concerted 
action of two bases - Tyrl53 and Lys65. Protons are trans- 
ferred by the proton relay system from the nucleophile to 
Tyrl53 and from Tyrl53 to Lys65, and a nucleophilic attack 
occurs, followed by the formation of the second tetrahedral 
intermediate. 

At the decomposition of the second tetrahedral intermedi- 
ate followed by the release of the second reaction product, 
the Lys65 residue cedes a proton to the Ser62 oxyanion and 
the enzyme returns to its initial state. 
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CONCLUSIONS 

The bioinf ormatic and phylogenetic analysis of the penicillin- 
binding protein family including D-aminopeptidases was car- 
ried out, and the conserved residues were identified. Three of 
them - catalytic Ser62, Lys65, and Tyrl53 - form a catalytic 
triangle - a specific proton relay system that captures/cedes 
a proton in the course of catalytic events and makes possible 
the activation of the generalized nucleophile in the D-amin- 
opeptidase catalysis. Molecular modeling showed that a con- 
served residue (Lys65) possess an unusually low pKa value, 



which was confirmed by an experimental pH-profile of the 
D-aminopeptidase catalytic activity. The specific role of the 
Lys65 residue in a catalytic mechanism of D-aminopeptidase 
was proposed at the formation and hydrolysis of the acyl- 
enzyme intermediate. • 

This work was supported by the Russian Foundation for 
Basic Research and Japanese Society for the Promotion of 
Science (grant M 09-08-92104). 
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